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ABSTRACT 



We investigate neutrino processes for conditions reached in simulations of 
core-collapse supernovae. Where neutrino-matter interactions play an important 
role, matter is partially degenerate, and we extend earlier work that addressed 
the degenerate regime. We derive expressions for the spin structure factor in 
neutron matter, which is a key quantity required for evaluating rates of neutrino 
processes. We show that, for essentially all conditions encountered in the post- 
bounce phase of core-collapse supernovae, it is a very good approximation to 
calculate the spin relaxation rates in the nondegenerate limit. Wc calculate spin 
relaxation rates based on chiral effective field theory interactions and find that 
they are typically a factor of two smaller than those obtained using the standard 
one-pion-exchange interaction alone. 

Subject headings: dense matter; neutrinos; stars: massive; supernovae: general 
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Introduction 



Understanding tlie mechanisms responsible for supernovae and neutron star formation 
requires a knowledge of the equation of state and transport mechanisms in matter at 
densities of the order of that in atomic nuclei and at temperatures ranging up to lO^^K. 
Despite almost half a century of work on the subject, the question of how a fraction of the 
large thermal energy from the collapse of the core is transferred to the outer stellar layers, 
thereby causing a supernova explosion, is one that has yet to find a convincing answer. 
Because of the high matter density in the new-born protoneutron star, most macroscopic 
transport processes are ineffective, and a variety of other mechanisms have been considered. 
These include energy transfer by neutrinos that interact with matter via weak interactions 
( [Colgate fc White|[l966l [Bethe fc Wilsol^[l985| ) , magnetic fields in combination with rotation 



(Ardeljan et al. 2004 Kotake et al. 2006), convection (Herant et al. 1994), or pressure 



waves and Alfven waves that may be emitted by the protoneutron star (Burrows et al. 



2007 Sagert et al. 2009). Some of these mechanisms lead to very characteristic features 



in the gravitational wave (Ott 2009) and neutrino signatures (Dasgupta et al. 2010), and 



these will be constrained by observations of the next galactic supernova. Other mechanisms 
have been tested in simulations, which demonstrate that, with current microphysical input. 



neutrino transport alone is insufficient to generate explosions ( Liebendorfer et al. 2001 



Rampp et al. 2002 Thompson et al. 2003). However, it is possible to obtain explosions 



when neutrino heating and fiuid instabilities are combined in axisymmetric models of stellar 



core collapse (Buras et al. 2006; Marek & Janka 2009 Bruenn et al. 2009 Suwa et al. 2010 



Brandt et al. 2011). 



In all the above mechanisms, transport by neutrinos is crucial: it is the dominant 
process for energy emission, it contributes to the transport of energy and lepton number, 
it infiuences the radial entropy and lepton fraction gradients that determine the stellar 
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structure and fluid instabilities, and it gives rise to the neutrino luminosities and spectra, 
which will become the main observables for probing the properties of matter at high 
density in the next nearby supernova event. Calculations of rates of neutrino processes 
in dense matter shortly after the discovery of weak neutral currents were reviewed by 



Freedman et al. (1977). In the early work, the particles participating in the reactions were 



taken to be free, but, subsequently, effects of strong interactions were taken into account 



(Sawyer 1975; Iwamoto & Pethick 1982). More recently, detailed calculations of rates of 



neutrino processes have been performed within a mean-field approach (the random phase 



approximation) by Reddy et al. (1998); Burrows & Sawyer (1998); Reddy et al. (1999); 



Burrows et al. (2006). One effect not included in these calculations is that excitations can 



decay due to interactions in the dense medium. As stressed in Raffelt & Seckel (1995); 



Raffelt et al. (1996); Hannestad & Raffelt (1998), this can lead to an energy transfer in 



neutrino processes considerably greater than that predicted on the basis of excitations with 



infinitely long lifetimes (i.e., from recoil effects alone). Lykasov et al. (2005) showed how to 



include these damping effects in a mean-field approach, and a unified approach to structure 



factors was described by Lykasov et al. (2008). Detailed calculations were performed based 



on chiral effective field theory (EFT) interactions in Bacca et al. (2009) for degenerate 
neutrons. The prime purpose of the present paper is to extend these studies to partially 
degenerate and nondegenerate conditions. 

The plan of the paper is as follows. In Section |2| we analyze the results of simulations 
of core-collapse supernovae and determine the conditions for which it is important to 
know rates of neutrino processes. These results point clearly to the need for a better 
understanding of neutrino properties in regions where nucleons are partially degenerate or 
nondegenerate. In Section |3| we develop a general formalism based on Landau's theory 
of normal Fermi liquids for calculating structure factors of strongly interacting matter. 
The spin relaxation rate for partially degenerate conditions is derived in Section |4] and the 
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nondegenerate limit is studied in Section 4.1 A key quantity is the spin relaxation rate 
in partially degenerate neutron matter, and we calculate this in Section [5] based on the 
one-pion-exchange approximation for nucleon-nucleon interactions, which is the standard 



one used in simulations (Hannestad & RafFelt 1998), and from chiral EFT interactions. 
Particular attention is paid to conditions of importance in supernova simulations. Finally, 
we summarize and give future perspectives in Section |6l 



2. Physical conditions in stellar collapse 

We consider the results of numerical simulations of the collapse of 15Mq and AOMq 
stars in a spherically symmetric supernova model with Boltzmann neutrino transport 
(Fischer et al. 2009). Figure [l] shows an overview of the conditions that are reached after 



the collapse of the stellar core of a ISMq progenitor star. The upper panels are constructed 
from a two-dimensional histogram in the parameter space spanned by the matter density p 
and the matter temperature T. The color of each bin [pj,pj+i], [Tj, Tj+i] is chosen according 
to the value of a time- and mass-weighted measure of the occurrence of the conditions 
corresponding to the bin, where the labels on the color bar specify the value of 



r-M 

if(p,r) = iogio' 



7^ dt dMe„clXij(i^,Mencl) , (1) 

^^0 Jo Jo 



where Xijity Mend) = 1 if the density and temperature are reached in the bin. 
Pi ^ pit, Mend) < Pi+i and Tj ^ T{t,Mend) < ^i+i, and Xij(t, M^nd) = otherwise. The 
integral over time t runs from zero (at the time of core bounce) to the final time r, while 
the integration over the enclosed mass Menci(-R) = dr Aur"^ p{r) runs from zero at the 
center of the star to the total mass M considered in the computational domain. Hence, a 
white bin in the p-T-plane indicates that the corresponding conditions rarely occur in the 
simulation, while a red bin indicates conditions occurring for a large mass and/or for a long 
time duration. The lower panels are constructed in the same way, but with the temperature 
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in Equation ([T]) replaced by the electron fraction Y^. The symbols indicate where neutrinos 
decouple from matter. This was determined by comparing the mean neutrino energy in 
the simulation with the mean neutrino energy expected for neutrinos in equilibrium at the 
temperature of the matter. The symbols mark the conditions under which the difference 
between the two mean energies was equal to 



70. 



Figure [2] shows the conditions for a simulation that was launched from a more massive 
progenitor star with a main sequence mass of 40Mq. All panels show on the far right-hand 
side the cold core of the protoneutron star at densities around 3 x lO^^gcm^^, temperatures 
around 10 — 20MeV and an electron fraction just below = 0.3. The black solid line 
shows the baryon Fermi temperature as a function of density (see caption) and thus the 
neutrons in this regime are degenerate. This matter lies in the inner part of the stellar core, 
which is collapsing subsonically, and it is unaffected by the bounce shock that forms outside 
the core at the sonic point and runs outward. The compactness of this highest density 
regime is sensitive to the equation of state (EOS) of neutron-rich matter. In parametric 
EOSs, this depends mainly on the incompressibility and the symmetry energy. If the EOS 
is stiff, the mantle of the protoneutron star will lie at a larger radius in a lower gravitational 
potential. This leads to comparatively small neutrino luminosities and soft spectra. If the 
EOS is soft, the core of the protoneutron star will be more compact and the mantle will 
lie at a smaller radius deep in the gravitational potential. This leads to larger neutrino 
luminosities and harder spectra. This general relationship between neutrino emission and 
proto-neutron star size applies to all factors that influence the compactness of the core: 
general relativistic effects ( Liebendorfer et al. 2001 Bruenn et al.||2001 ), different progenitor 



masses (Liebendorfer et al. 2003), and different EOSs (Sumiyoshi et al. 2007). In the left 



panels of Figures [T] and |2| we show results from simulations based on the Lattimer-Swesty 
EOS ( Lattimer fc Swesty| 1991), which produces a more compact protoneutron star core. 



while the right panels show results from a simulation with the stiffer Shen et al. (1998) EOS. 
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New constraints from neutron matter calculations (Hebeler et al. 2010) and from modeling 



X-ray burst sources ( Steiner et al.||2010 ) favor EOSs that yield more compact (cold) neutron 



stars than the Shen et al. (1998) EOS can give. 



When the bounce shock runs through the mantle of the protoneutron star, it causes a 
temperature jump. Hence, in the upper panels of Figures [T] and [2} matter in the cold core 
is well separated from that in the shock-heated mantle, which is found above the black line 
that indicates partially degenerate conditions. The conditions in the mantle are represented 
by a rather narrow band in the p-T-plane that runs from densities around lO^^gcm"^ at 
temperatures of ~ 30MeV down to densities of lO^^gcm"^ at temperatures of ~ 3MeV. 
This band is approximately described by the relation 

/ \ 1/3 

T = 3 ( ——^ r I MeV , (2) 

VlOiigcm-^ y ' ^ ' 

for which the temperature at p = lO^^gcm"^ is equal to the Fermi temperature at 
that density. In Figures [T] and [2| this relation is shown as dashed lines. This behavior 
may be understood as being a consequence of the fact that the shocked matter expands 
approximately adiabatically (neglecting neutrino heating). Due to the high temperature, 
most of the entropy resides in the leptons and the photons, which behave as ideal, relativistic 
gases, for which the entropy per unit volume varies as T^. Since the baryon density is 
proportional to p, the temperature on adiabats behaves approximately as T ~ p^/^. Note 
that the maximum temperature depends strongly on the mass of the progenitor star and, 
in the case of the ISM© model, also on the EOS. During the first second of the post-bounce 
phase, more than 70% of the z/g, Ve and v^^r luminosities are emitted from this hot mantle of 
the protoneutron star. At higher densities, the influence of the neutrinos on the supernova 
dynamics is comparatively small, because neutrinos are trapped so that thermal and weak 
equilibrium with the matter is established. At very low densities, the influence of neutrino 
processes is limited by small weak-interaction rates. The sizable electron chemical potential 
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in the mantle favors electron capture over positron capture, so that it is the neutrino 
emitting protoneutron star mantle that deleptonizes most efficiently and becomes the most 
neutron-rich region in the computational domain. The resulting low electron fractions 
in the density interval p = lO^^gcm"^ to lO^^gcm"'^ can be seen in the lower panels of 
Figures [T] and [2] 

Neutrino reaction rates exhibit a strong energy dependence; as a consequence, neutrinos 
with different energies decouple from the matter at different locations. As an approximate 
measure of where decoupling occurs, we take the radius at which the mean neutrino energy 
obtained from the simulation begins to differ from the one expected for a gas of neutrinos 
in thermal equilibrium with the matter. The symbols in Figures [T] and |2] indicate the 
thermodynamic conditions at these points. To indicate the time evolution, the conditions for 
decoupling are shown at 50 ms after bounce (open symbols) and at the end of the simulation 
(filled symbols). In general, the neutrinospheres for electron neutrinos (square symbols) are 
found at the lowest density due to their charged-current reactions with neutrons, which are 
the most abundant species under the given conditions. Electron antineutrinos, which have 
charged-current reactions with protons in addition to neutral-current reactions, decouple on 
average at densities ten times larger (diamonds). The /i and r neutrinos and antineutrinos 
interact with matter only via neutral-current reactions and therefore decouple at the highest 
density, as indicated by the stars. This density hierarchy translates into a temperature 
hierarchy for the regions from which neutrinos are emitted. As a result, the spectra of r 
and u^^r are harder than that of Ug, which in turn is harder than that of Uf.. 

The post-bounce evolution can be separated into two phases. During the ffist 
phase the mass of the protoneutron star is significantly lower than the maximum mass 
supported by the EOS. Thus, the protoneutron star does not contract significantly and the 
neutrinospheres move to higher density with time as the accretion rate and the density of 



- 9 - 



the accreting layers decrease. Because regions at higher density have higher temperature, 
this also leads to a continuous increase of the temperature of emitted neutrinos. However, 
in spherically symmetric models it can happen that the steady cooling by z/^,,- and Uf^^r 
emission within a narrow range of the mass coordinate decreases the local temperature 
sufficiently that the mean energy of z/^^^ and Uf^^r decreases and eventually falls below that 



of the z/e (Liebendorfer et al. 2004). The first phase of the post-bounce evolution is well 
represented by the ISM© model shown in Figure [1} After the onset of the supernova 
explosion, which was not successful for the models shown in Figures [T] and |2j the mean 



energies of the neutrino spectra decrease as the protoneutron star cools ( Fischer et al.|[2009 ). 



Due to the lower opacity for neutrinos with lower energy, the neutrinospheres continue 
to shift to higher density and make the high-density part of the protoneutron star more 
' neutrino- visible ' . 

The first phase of the postbounce evolution dictates the success of the supernova 
explosion, i.e., whether sufficient entropy can be accumulated behind the standing accretion 
shock to revive it and drive it through the outer layers of the star. This phase of the 
supernova mechanism is extremely challenging to model because the emission, transport 
and absorption of neutrinos behind the shock couples to vigorous three-dimensional fluid 
instabilities between the neutrinospheres and the accretion shock. These multidimensional 
effects are not included in the data shown in Figures [T] and [2] 

If an explosion has not taken place during the first stage, the post-bounce evolution 
enters a second phase when the mass of the protoneutron star approaches the maximum 
stable mass. This scenario may occur in the rarer cases of 'failed supernovae' that are 
triggered by the collapse of a very massive progenitor star. As the limiting stable mass 
of the protoneutron star is approached, accretion makes the protoneutron star more 
compact and pushes matter with high entropy into the regime of the neutrinospheres. The 
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corresponding temperature change increases the opacity so that the neutrinospheres shift to 
lower densities, as can be seen in Figure [2] This leads to a dramatic increase in the mean 



energy of emitted u^^r and Uf^^r (Liebendorfer et al. 2004 Sumiyoshi et al. 2007). However, 



only the AOMq model reaches this second phase within the duration of the simulation. With 
the Lattimer-Swesty EOS, it takes 0.43 s to form a black hole, while for the stiffer Shen et 



al. EOS it takes 1.4s (Sumiyoshi et al. 2007 Fischer et al. 2009). 



In summary, the emission of neutrinos from the partially degenerate hot mantle of 
the protoneutron star, the EOS at nuclear densities, and the accretion rate set by the 
progress of the supernova shock and the density of the outer layers of the progenitor star 
are the key ingredients that determine the compactification of the protoneutron star during 
the early post-bounce evolution. The emission of neutrinos also shapes the radial entropy 
and lepton number profiles that determine where stratified layers of matter are stable 
against convection. The emitted neutrinos finally can be absorbed by matter at larger 
radii between the protoneutron star mantle and the standing accretion shock. The heat 
transferred per absorbed neutrino increases with the difference between the temperature of 
the neutrinos and that of the absorbing matter. This difference increases with the distance 
from the neutrinospheres, but the neutrino number flux decreases with the distance because 
of geometrical dilution. Hence, there is a region close to the protoneutron star where 
heating is most effective. This peak heating at the base of the matter lying between the 
protoneutron star and the standing accretion shock establishes a negative entropy gradient. 
The entropy gradient induces fluid instabilities that start to transport energy to the shock 



by mechanical turnover (Herant et al. 1994 Marek & Janka 2009). The perturbations in 



the convective fluid flow are enhanced by the standing accretion shock instability (Blondin 



et al. 2003 Foglizzo et al. 2007). The fluid instabilities and the net rate of neutrino cooling 



and heating determine the further evolution and geometry of the supernova. It is therefore 
essential for supernova models to accurately treat the neutrino-matter interactions in the 
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hot mantle of the protoneutron star under the conditions shown in Figures [T] and |2j which 
we shall investigate in Section [5] 

3. General formalism 

Motivated by the relevant supernova conditions, we extend our work on neutrino 



processes in degenerate matter (Lykasov et al. 2008 Bacca et al. 2009) to the case of 



partially degenerate and nondegenerate neutrons. 

The basic ingredients in the calculation of rates of neutrino processes are the structure 



factors for the vector and axial vector response (Raffelt 1996). For nonrelativistic neutrons, 
the vector current couples to the neutron density and the axial current to the neutron 
spin. We denote the coupling strengths by Cy and Ca, respectively, and for free neutrons 
Cy = 1/2 and Ca = fi'A/2 = 1.26/2. We shall focus on the axial current response for two 
reasons. The first is that it is responsible for the largest contributions to rates of neutrino 
processes. The second is that axial current processes have potential for equilibrating 
neutrino energy distributions because the axial current is not conserved and therefore it is 
possible for the energy transfer to or from the neutrinos to be nonzero even for processes in 
which there is a small momentum transfer. For the vector interaction, the energy transfer 
vanishes for zero momentum transfer because the vector current is conserved. For neutrons, 
the structure factor for the axial current therefore has the form 

Sx{uj,ci) = C\S^{uj,(^) , (3) 

where the spin structure factor is given by 

S„{uj,(\) = —- ^-^ImXa(w,q) • (4) 

Here Xo- is the spin-density-spin-density response function. (We use units in which 

h = c = k-Q = 1.) 
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We begin by presenting a phenomenological approach. As has been stressed by 
Raffelt and coworkers, it is a good first approximation in many situations to consider the 
long-wavelength response because the momentum transfer to the nucleons is of the order of 
a typical neutrino momentum, which is generally small compared with the characteristic 
momentum of a nucleon. The latter is of the order of the Fermi momentum (Svr^n)^/'^ for 
degenerate matter and the thermal momentum (mT)^/^, where m is the nucleon mass, for 
nondegenerate matter. We consider the relaxation of S, the total z-component of the spin 
of the system to its equilibrium value S^q, which is given in linear response theory by 

Scq = —XM Uz , (5) 

where —11^ corresponds to the strength of a magnetic field that polarizes the system and xm 
is, apart from factors, the static magnetic susceptibility. If one assumes that the approach 
of S to its equilibrium value is proportional to the difference between S and its equilibrium 
value, we may write 

dS S — Sen ,„s 

-Tt= 

where r is a characteristic relaxation time. The quantity 1/r corresponds to what 



Hannestad & Raffelt (1998) refer to as the spin fluctuation rate. The total spin of the 
system thus approaches its equilibrium value exponentially, if the magnetic fleld does not 
depend on time. 

On Fourier transforming Equations ([s]) and ^ one flnds for the frequency-dependent 
response function 

^-^"^ = -Z4 = r^- 

In order to make detailed calculations, we shall consider nucleonic matter as a system of 
interacting quasiparticles. We shall assume that the spectrum of a single quasiparticle may 
be described by an effective mass m* that is independent of its momentum p, and therefore 
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the energy of a single quasiparticle is given by 

6p. = ^o + £^. (8) 

We take into account two-body quasiparticle interactions, with a spin-dependent part of 
the form 

/po-i,p'cr2 = 5'o CTi ■ cr2 , (9) 

where cTj are Pauli matrices. The theory may be regarded as an extension of Landau's theory 
of Fermi liquids generalized to higher temperatures when the quasiparticle interaction is 
independent of the momenta of the quasiparticles. The parameters m* and Qq generally 
depend on temperature and density, and their values are chosen to ensure that bulk 
properties such as the specific heat and magnetic susceptibility agree with the results of 
more microscopic calculations. 

We shall take the coupling of the spin to the external field Uz to be given by a change 
in the energy of a quasiparticle, ep^, equal to 

5ep^ = aUz. (10) 



By generalizing the standard calculation of the static magnetic susceptibility (Baym & 



Pethick 1991 ) to higher temperatures, one finds 





1 + gox" 



where 

_ dn{e 



r = -ET^. (12) 

per 



is the static susceptibihty in the absence of interactions between quasiparticles and 

= e(^-M)/T + 1 ' ^^^^ 
is the Fermi function, /i being the chemical potential. At temperatures low compared with 
the Fermi degeneracy temperature Tp = pp/2m*, where pp is the Fermi momentum, the 
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density of states x° = rn*pp/TT'^ for a single species of particles with two spin states, while 
at temperatures high compared with Tp, one finds the Curie law = n/T. On inserting 
Equation (12) into Equation ([T]) one finds 

.0 



X 



1 + 90X^ - ^^Ta 



where 



= r (1 + gox 



(14) 



(15) 



This equation has the same form as in Lykasov et al. (2008) [see, e.g.. Equation (24)]. The 
significance of the time To- may be brought out by writing Equation ^ in terms of the 
deviation of the spin from its local equilibrium (I.e.) value 



(16) 



which is the equilibrium spin density in a system exposed to an external field ~Uz plus 
a "molecular field" go S due to interactions with the other particles. Equation ^ then 
becomes 

S — Sie. 



dS 

~dt 



Trr 



(17) 



From kinetic theory, one saw in the degenerate limit how mean-field effects enter relaxation 



rates (Lykasov et al. 2008): the time does not depend on mean-field effects, while the 



characteristic time for relaxation of a deviation of the total spin from its equilibrium value 
is decreased by the factor by which mean-field effects decrease the magnetic susceptibility. 
The reason for this effect is that when the magnetic susceptibility is reduced, the free energy 
difference driving the return to equilibrium, which for a given spin deviation is inversely 
proportional to the magnetic susceptibility, is increased, and relaxation is faster. 



In the work of Hannestad & Raffelt (1998), the normalization of the spin response 



function was determined by the condition that it give correctly the static structure factor, 
which was assumed to be that of a noninteracting gas. It is therefore of interest to 
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investigate the relationship between the static structure factor and the spin response 
function when interactions are taken into account. The static structure factor 5*0- (q) is 
defined by 

/oo 
dcoS^{cu,q). (18) 
-oo 

Since the dynamical structure factor S{uj, q) is related to the response function by 
Equation (|4]), one has 

where we have used the fact that \m.Xa{'^i q) is an odd function of u. From the inequality 
|coth[w/(2r)]| ^ \2T/uj\, it follows that 

S,(q)>- r-5^^^^^^ = -X«, (20) 

n TT a; n 



the equality applying only if all the weight of Imx^. is at zero frequency. Using Equation ( 14 ) 
for the response function, one finds 

SM ^ 0) ^ - . (21) 
n 1 + ^oX° 



These results represent a generalization of Equation (15) of Hannestad & Raffelt (1998) 
to allow for spin correlations. For neutron matter at low temperature, QqX^ is equal to 
the Landau parameter Gq, which is calculated to be approximately 0.45 at a density 



of po/100 and 0.8 at nuclear saturation density, po = 2.8 x 10^^ gem ^ (Schwenk et al. 



2003), and therefore interactions depress the static structure factor by up to almost a 



factor of two. For density correlations, the equality sign applies in the long-wavelength 
limit, since total particle number is conserved and, consequently, matrix elements of the 



density operator to states with nonzero excitation energy vanish in this limit (Baym & 



Pethick 1991, Section 1.3.3). If interactions that do not conserve spin are neglected, the 



inequality (20) becomes an equality at long wavelengths. This approximation has been 
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employed in calculations of neutrino processes in nondegenerate matter at low densities 



based on the virial expansion (Horowitz & Schwenk 2006). 



A number of cautionary remarks about our approach should be made. First, in the 
above we have calculated only the contribution to the susceptibility due to the coupling 
of the external field to single quasiparticle-quasihole intermediate states. The effects of 
multipair states are taken into account in so far as they may be regarded as arising from 
single-pair states as a consequence of the collision term in the kinetic equation. Other 
processes that can result in the creation of two- and higher-quasiparticle-quasihole states, 



e.g., two-body contributions to effective operators (Menendez et al. 2011), have been 



neglected. As emphasized in Olsson et al. (2004), because of the noncentral character 



of nuclear forces, these are nonzero even for very long wavelengths of the exciting field. 
Similarly, we have neglected many-body contributions that can renormalize the magnetic 
moment (or in the case of the coupling to the weak field, the weak charge) of the nucleon. 
In estimating rates of processes one must multiply the above results by the square of the 
appropriate renormalization factors. 

One may ask how good is the assumption of a single relaxation time. This question 



was addressed for degenerate nucleons in Pethick & Schwenk (2009). There it was shown 



for u <^ T that, if the relaxation time is chosen to reproduce the leading behavior for 



t<j ^ I/to-, the real and imaginary parts of the response function given by Equation (14) 



differ from the result obtained by solving the transport equation exactly by less than 10%. 
For nondegenerate matter this question has not yet been addressed, but one can draw 
on experience with other transport problems for this case. The relaxation time at high 
values of lot^j is given by the simplest variational estimate for the relaxation time, which 
uses a trial function proportional to the quasiparticle spin. The actual relaxation time 
at low frequencies is longer than this estimate, but to date no estimates of it have been 
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made. For nondegenerate, spinless atomic gases one finds that relaxation times for viscosity 
and tfiermal conduction in tlie fiydrodynamic regime differ from tlie simplest variational 



estimates by amounts tliat are typically of order two per cent or less (Chapman & Cowling 



1970). The part of the interaction relevant for spin relaxation depends on the spin of the 



particles, and it has a different spatial dependence from typical interactions between spinless 
atoms, so it is desirable to calculate the ratio of the relaxation times in the hydrodynamic 
and coUisionless regimes for nuclear interactions. However, it seems unlikely that the 
deviations will be more than a few per cent, which is small compared with the uncertainties 
in the strong scattering amplitudes that enter the relaxation rates. For practical purposes 



it is therefore expected to be a very good approximation to use Equation (14) for all ut^. 



4. Relaxation time 

In this section, we calculate the response function microscopically for frequencies high 
compared to the collision ratej^ In terms of the phenomenological model discussed in the 
previous section, this is given by 

xM~i^. (22) 
Quantitatively, the response function may be calculated using field-theoretical methods. 



as described in Lykasov et al. (2005 ). The contribution to the response function proportional 
to l/cj follows from the diagrams shown in Figure |3] In the figure, the lines represent 
quasiparticles, and the interaction vertices are those for quasiparticles. The square boxes 
denote scattering vertices for quasiparticles, symmetrized with respect to the incoming and 



^Here we are referring only to the quasiparticle contribution to the response. By large 
frequencies, we mean frequencies large compared with l/r^, but still small compared to the 
other energy scales in the system. 
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outgoing states. The diagrams in Figure [3] (a) correspond to self-energy corrections which 
take into account the nonzero width of the quasiparticles, while those in Figure [s] (b)-(d) 
correspond to vertex corrections. 

As we have argued above, the leading contribution to the response functions at large u 
is imaginary, and a convenient way of calculating this is by use of the techniques developed 
in field theory ( Landau||1959 Cutkosky||1960 Taylor||1960 ), which have been exploited in the 



context of condensed matter in Langer (1961, 1962) and Carneiro & Pethick (1975). The 



energy denominators that lead to singularities may be identified by finding intermediate 
states in the response function which have the property that when lines corresponding 
to the intermediate state are cut, the diagram falls into just two parts. In Figure |3] such 
intermediate states corresponding to two-quasiparticle-two-quasihole states are denoted 
by dashed lines. By including all such intermediate states one automatically includes 
contributions that correspond both to "in-scattering" and "out-scattering" terms in an 
approach based on a kinetic equation. On evaluating the response function one finds 
27r 1 



= i ^ o y]^i(^i + - 0-3 - 0-4) 



1234 



X 



|(34|r|12)|V^2(l -ri3)(l -^4) - I (12|r|34)|'(l -ni)(l -712)713^4) 



(23) 



Here the index 1 is shorthand for picri and so on, ni are equilibrium distribution functions, 
5p = 5(pi + P2 — P3 — P4) and = 6{±uj + ei + 62 — 63 — €4), and we have dropped 
the momentum transfer q from the momentum conservation condition according to the 
arguments given above. The factor of 1/2 is to avoid double counting of final states 
when the sums over ps and p4 are unrestricted. The two different scattering matrices in 



Equation (23) may be combined if one makes use of the fact that the interaction is invariant 



under time reversal, which implies 



(34|r|12) 



-l,-2|r|-3,-4)^ 



(24) 
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where the minus sign indicates that both the momentum of a quasiparticle and its spin 
must be reversed. Since the distribution function depends on momentum only through its 
magnitude, since all momenta are summed over, and since each term is quadratic in the 



spins of quasiparticles, it follows that x{^) has the form of Equation (22) with 

- = ^ 5^ai(ai-a2-a3-a4)|(34|r|12)|'nin2(l-n3)(l-n4)5p [5+ - C] • (25) 

^'^ ^ 1234 

This form demonstrates explicitly the fact that the relaxation time is an even function of 
u, as is required by the condition that the real part of x be an even function of u and the 
imaginary part an odd function. By making use of the relation 

e-(^'-^)/^, (26) 



1 - rii 



and exchanging the roles of 1 and 2 with 3 and 4 in the sum in Equation (25), one can see 
that the second term is e~'^^^ times the first one, and therefore 

1 (l-e-'^/^) TT 



^ 1234 



J2 + ^2 - - a4)|(34|r|12)|' nin2(l - n3)(l - n,) 5p 5+ . (27) 



This may be expressed in the more symmetrical form 

1 2 sinh[(X'/(2T)l vr v-^ / x i - i^i, i2 ,- . -i-r 1 

(28) 

To convert the sums over momenta to integrals, we use as momentum variables the total 
momentum P = pi + p2 = P3 + P4 and the relative momenta p = (pi — P2)/2 and 
p' = (p3 — P4)/2. Because all directions of the momenta are being integrated over, it is 
convenient to work with the quantity 

W = ^ E Tr[(34|r|12ra^[(^l + ^2)^(34|r|12)]], (29) 

fc=l,2,3 

which depends only on P,p, and p', and three angles specifying the relative orientation of 
the vectors P, p, and p'. The trace Tr is over the spin indices for all the quasiparticles and 
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the spin structure with the commutator is the same as for degenerate conditions (Lykasov 



et al. 2008). The integration over may be ehminated by making use of the energy 



conservation condition, 



The final resuh is 



^-='-^u. (30) 



1 1 2sinh[u;/(2T)] m* /"^ ''^^ 
— = — T / a cos I a cos 

2^6 ^ ^0 ' 



1 J -I JO 

4 



POO POO ^ -I 

where 6' is the angle between P and p, 6' the angle between P and p', and is the angle 
between the normal to the plane containing P and p and the normal to the plane containing 
P and p'. For definiteness, we have taken u to be positive. For negative u, the lower 
limit of the integral over p is {m*\uj\Y^'^ . The integral is symmetrical in p and p' since 
2dpp^ ijp' + m*ujy^'^ = pp'dp^ = pp'dp''^. 

Alternative forms for the relaxation rate that bring out more clearly the physical origin 
of the various contributions are obtained by replacing the factors 



in Equation (31) by 

(p2 + m*ujy/^ W [n,n2{l - n3)(l - n^) - (1 - ni)(l - n2)n3n4] , (33) 

or by 

[(p^ + m*ujy/^ - (p2 - m*ujy/^] Wnin2{l - n3){l - n^) . (34) 
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4.1. Nondegenerate limit 



In the nondegenerate limit, /x/T — t- — oo, the product of thermal factors in Equation (31 ) 
becomes 

which shows that the distribution functions for the relative momentum and the center-of- 



mass momentum are independent of each other. However, the integral in Equation (31) is 
not particularly simple, because W is generally a function of all 5 integration variables. 
Simplification is possible if the transition probability W is independent of the center-of-mass 
momentum, which is the case if the influence of the medium on the scattering process is 
neglected. Then is a function of p, p' = {p^ + m*u;)^/^, and 0, the angle between p and 
p'. The expression for the relaxation rate then takes the form 

8ir5 uj x" 7-1 Jo 

POO 

Jo 

An uj/T 



where 



{■■■) = — —^ , (38) 



2 



/'OO n 1 

j cipp^e-P'/^""'^^ J dcosG.... (39) 



/7r(m*T)3/2 Jg 

In the nondegenerate limit, the dependence on density factors from the integral. We 
therefore introduce the quantity 

^ 2sinh[a;/(2T)] _^ ^ ^-./i2T)^ ^ (4q) 

u/T 
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in terms of which the spin relaxation rate is given by 



1 m*n ^ pm*/m„ 



5. Results 

We next calculate the spin relaxation rates using nucleon-nucleon (NN) interactions 



based on chiral EFT to next-to- next-to-next-to-leading order (N^LO), as in Bacca et 



al. (2009) for degenerate conditions, and compare with the results obtained from the 



one-pion-exchange (OPE) approximation to nuclear interactions, which provides the 



standard rates for bremsstrahlung in supernova simulations (Hannestad & Raffelt 1998). 
At this level, the transition amplitude is independent of the center-of-mass momentum 
and given by the antisymmetrized interaction (12|T|34) = (12|(1 — Pi2)^n|34) (with 
exchange operator Pu). To be explicit regarding the normalization, the direct and exchange 
contributions for the OPE approximation are given by 

^ ^<Ti ■ ko-2 • k <Ti ■ k' 0-2 ■ k' + fc'2(l - (Ti ■ o-2)/2" 



;12|(1-Pi2)1^nT|34) = -(^ 



f^2 _|_ ^2 _|_ 



(42) 



with pion decay constant = 92.4 MeV, neutral pion mass m,r = 134.98 MeV, and 
momentum transfers k = pi — p3 and k' = pi — p4. For the calculation based 
on chiral EFT interactions, we make a partial-wave expansion with the convention 
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p) = X]im^^^^/m(p) After coupling to spin, this leads to 



2sinh[a;/(2T)] A (An) 



C 



[ml — msm'g) 



^{m*T)3/^ u/T 

i i' L]j i' i L 
1 1 J J 1 1 1 J 

POO 

X / dpp" + rn*a;e-P'/(™*^)-"/(2^) {^/p^ + m*u \V^T^\p){^/p^ + m*u |V//='|p) , 

(43) 



i L 
V 



with standard notation for the quantum numbers as in Lykasov et al. (2008), and where the 



sum is over allowed partial waves with matrix elements ^p\yli\\pl). 

We begin by comparing results in the OPE approximation allowing for partial 
degeneracy with those for the nondegenerate and degenerate limits. For all results in the 
following, the neutron effective mass m* is taken to be the bare mass m. In Figure [4], we 
show the spin relaxation rate as a function of temperature for densities p = lO^^gcm"^ (top 
panel) and p = lO^'^gcm"^ (bottom panel). The full (partially degenerate) results are based 



on numerical integrations of Equation (31). The spin relaxation rate in the degenerate 



limit (Lykasov et al. 2008) is shown for T ^ Tp/n and the nondegenerate limit is given 



by Equation (37). Our results clearly demonstrate that the nondegenerate limit is a good 



approximation to the full results down to temperatures well below the Fermi temperature, 
Tp = 30MeV at a density p = lO^^gcm'^, and Tp = 6.5 MeV at p = lO^^gcm"^. In the 
nondegenerate case, the spin relaxation rate is linear in the density, and information about 
the temperature dependence is contained in the quantity r.{uj), which is plotted in Figure |5 



for cu = 0. As expected from the spin relaxation rate in the degenerate limit (Bacca et al. 



2009), results based on chiral NN interactions at N^LO are typically a factor of two smaller 



than those obtained using the standard OPE approximation. As in Bacca et al. (2009), we 



consider different N^LO NN potentials (Entem & Machleidt 2003 Epelbaum et al. 2005) 



and the resulting band gives a range of uncertainty at this level. We also give a fit to the 
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central value of the N^LO band: 



0) = 0.133 i-f- 



1.12 



-0.02(T/MeV) 



(44) 



Since at typical frequencies and temperatures, the oj dependence is weak (see Fig. [T]), we 



recommend substituting the w = fit, Equation (44), to calculate the spin relaxation rate. 



Equation (41). This can be combined with Equations (|4]) and ([T]) for the spin dynamical 
structure factor to explore our rates for neutrons in astrophysical simulations. 

In addition, we compare the N^LO results to those based on the T matrix for NN 



scattering, which has been used for degenerate neutrons in Hanhart et al. (2001). For 
w = 0, the partial- wave potential matrix elements {p\V^l,~^\p){p\V^l,~'^\p^ in Equation (43) 



are replaced by on-shell T matrices, {p\Ti^, ^\pY (j)\Ti^i (note the complex conjugate). 



which are give in terms of scattering phase shifts and mixing angles [see, e.g.. Brown 



& Jackson (1976)]. For NN scattering, we take the phase shifts and mixing angles 



from NN-OnLine (2011) for ii^iab ^ 350 MeV, and note that contributions from higher 
i?iab = 2]9^/m can become important for temperatures T > 30 MeV. For temperatures 
T < 30 MeV, the T matrix results are similar, but somewhat lower compared to those from 
chiral NN potentials at N'^LO. This indicates that noncentral neutron-neutron interactions 
may be perturbative in chiral EFT for the energies relevant to the spin relaxation rate 
(where low energies are also suppressed by one power in momentum p at w = 0). 

We now consider how good the expression for the spin relaxation rate in the 
nondegenerate limit is under relevant conditions encountered in simulations of core-collapse 
supernovae. As representative of conditions during outflow, based on Figures [T] and [2} we 
take for values of the density ranging between po/100 and po the temperature given by 
Equation (g. 



In Figure |6| the full (partially degenerate) results are compared to the nondegenerate 
limit for these conditions. This figure demonstrates that results for the nondegenerate limit 
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are an excellent approximation for matter in supernova simulations at subnuclear densities. 
Moreover, for the broad density range considered, the results based on N'^LO NN potentials 
are similar to those based on the NN T matrix, and both lead to spin relaxation rates that 
are a factor two or more smaller compared to the OPE approximation. 

Next, we study the dependence of the spin relaxation rate in the nondegenerate 
limit on frequency. When u ^ 0, the matrix elements that enter the rate are off-shell 



for p ^ p' = a/p^ + m*uj [see Equation (43)], because the energies of the particles in the 
initial and final states differ by to. Figure [T] shows as a function of co/T for the OPE 
interaction for three different temperatures, and compares the exact u dependence with 
results in the on-shell approximation, where the partial-wave matrix elements are evaluated 
at a mean energy V^(a/p^ + 'm*u/2, ^/p^ _|_ rn*u]/2). The differences are found to be small, 
because the momentum dependence of the interaction matrix elements varies on a scale of 
the order of the pion mass (and higher momenta for chiral EFT interactions), which is large 
compared to the energy transfer, which is of the order of the temperature. This conclusion 
is reinforced by the plots of the relative difference of the spin relaxation rate calculated with 
the exact u dependence compared to the on-shell approximation in Figure [8] for typical 
conditions in supernova simulations. 

For astrophysical applications, the quantity that enters the rates of neutrino processes is 
the spin dynamical structure factor which is related to the relaxation rate by Equations (|4]) 



and (14). In Figure [9| we plot the structure factor for two different densities, 10^'^ gcm~^ and 
lO^^gcm^^ [with corresponding temperatures given by Equation ([2])]. The blue bands show 
the results for the N'^LO NN potentials and the solid lines those for the OPE approximation 
calculated within the formalism of this paper. As one would expect from the fact that 
the width of the peak in the structure factor is proportional to the relaxation rate and its 
height is inversely proportional to the rate, the results for the N^LO interactions have a 
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higher peak value and a narrower width compared with those for the OPE approximation. 
As representative of earher work, we show as dotted hnes results for the OPE calculations 



of Hannestad & Raffelt ( 1998 ) as implemented in the Basel code for core-collapse supernova 
simulations (Fischer et al. 2012). For a density of lO^^gcm"^, the Hannestad & Raffelt 



(1998) results lie 30-60% below our OPE results, while at 10^'^ gem ^, the difference is 



generally less. The difference between the OPE results may be a consequence of imposing a 



normalization condition in Hannestad & Raffelt (1998), while in our approach the width is 



obtained consistently by solving the Boltzmann equation. For the comparison in Figure [9} 
we have neglected the frequency dependence of the relaxation rate, which is weak at these 
temperatures and would increas the structure factor with increasing uj/T (see Figure [?]), 



as well as Fermi liquid corrections [g^ = in Equation (14)], which are not included in 



Hannestad k Raffelt (1998). 



In scattering, the energy transfer to or from neutrinos is comparable to the width 
of the peak in the structure factor, and thus one sees that the energy transfer under 
realistic supernova conditions can be comparable to the temperature. It is important in 
future simulations to allow for this energy transfer. In addition, the rate of neutrino pair 
bremsstralung is roughly proportional to the spin relaxation rate, and therefore the N^LO 
results predict a lower rate for this process than does the OPE approximation. This is 



consistent with what has previously been found for degenerate conditions (Bacca et al. 



2009). More comprehensive comparisons of predictions for rates of neutrino processes will 



be presented in future work that will extend the present calculations to take into account 
the effect of protons. 
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6. Concluding remarks 



We have analyzed simulations of core-collapse supernovae and have shown that in the 
post-bounce phase the relevant conditions for neutrino processes are partially degenerate or 
nondegenerate. We then developed a formalism for calculating neutrino rates in strongly 
interacting matter by generalizing to the partially degenerate regime the approach used for 



degenerate matter in Lykasov et al. (2008). The resulting spin dynamical structure factor 



takes into account both mean-field effects and collisions between excitations in neutron 
matter. We then calculated the spin relaxation rate, a key ingredient in the structure 
factor which enters expressions for the rates of neutrino processes. This was done at two 
levels of NN interactions, the OPE approximation, which is commonly used in simulations 



(Hannestad & Raffelt 1998), and chiral EFT interactions. We have found that chiral NN 



interactions lead to a reduction of the spin relaxation rate typically by a factor of two for a 
broad range of conditions. This reduction is similar to what previously has been found in 



the degenerate regime (Bacca et al. 2009). We have also found that our OPE rate, where 



the width is obtained consistently by solving the Boltzmann equation, differs from the OPE 



results of Hannestad & Raffelt (1998). This may be a consequence of the imposition of a 



normalization condition in Hannestad & Raffelt (1998). Moreover, our results demonstrate 



that the nondegenerate limit is an excellent approximation for the conditions encountered 
in the post-bounce phase of matter at subnuclear densities in supernova simulations. 

Future directions include the study of many-body contributions and of many-body 
forces and electroweak currents in chiral EFT, and the extension of the calculations to 
mixtures of neutrons and protons, where the central parts of nuclear interactions can cause 
relaxation of the axial charge, because of the different axial charges of the neutron and 
proton. These extensions will be greatly simplified by the finding that the nondegenerate 
limit provides an excellent approximation for the relevant supernova conditions. 
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Throughout the paper we have assumed that the basic couphng of the weak neutral 
field (that of the Z boson) to nucleons is via a one-body operator. However, the spin 
relaxation effect that we have considered amounts to coupling a single quasiparticle- 
quasihole pair to two quasiparticle-quasihole pairs. In other words, strong interactions 
have generated a two-body contribution to the weak charge operator. However, not all 
two-body contributions to the operator are included in calculating relaxation effects by 
this procedure. The approximation we have made includes only contributions with a single 
quasiparticle in an intermediate state, which are enhanced by Consideration of other 
two-body contributions to the weak charge operator is left for future studies. 
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Temperature: Cosmic Matter in the Laboratory", and by the DFG through Grant SFB 
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Fig. 1. — Color-coded plots of the occurrence of thermodynamic conditions during the post- 
bounce phase of spherically symmetric supernova simulations launched from a 15Mq pro- 
genitor star. Left (right) panels refer to a simulation with the Lattimer-Swesty (Shen et 
al.) equation of state (EOS). Upper (lower) panels show temperature T (electron frac- 
tion Ye) versus mass density p. The occurrence is expressed on a logarithmic scale by the 
amount of matter which experiences a specific set of conditions, normalized to Mq, inte- 
grated over time during the post-bounce phase [see text and Equation ([T]) for details]. The 
solid lines indicate where the temperature is equal to the Fermi temperature of the baryons, 
Tp = (37i'^p/m)'^^^ /2m, and the dashed lines represent the conditions specified by Equa- 
tion In each panel, six different symbols indicate the conditions under which neutrinos 
decouple thermally from matter. Squares refer to z/g, diamonds to z/g, and stars to u^ j. 
and Ufj^^r- Unfilled symbols refer to the early post-bounce phase 50ms after bounce), 
whereas filled symbols refer to the latest times considered (~ 400 ms after bounce for the 
Lattimer-Swesty EOS simulation and ~ 430ms after bounce for the Shen et al. case). 



- 30 - 




Fig. 2. — Same as Figure [T| but for a 40Mq progenitor. As in Figure [T| unfilled symbols 
refer to the early post-bounce phase 50ms after bounce), whereas filled symbols refer to 
the latest times considered, which is close to black-hole formation 430 ms after bounce 
for the Lattimer-Swesty EOS simulation and ~ 1.36s after bounce for the Shen et al. case). 
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Fig. 3. — Diagrams for the imaginary part of the response function. Sohd Unes represent 
quasiparticle propagators and the square boxes interactions of quasiparticles. The dashed 
hues represent two-quasiparticle-two-quasihole states that are on the energy shell. 
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Fig. 4. — Spin relaxation rate I/to- versus temperature T based on the one-pion-exchange 
(OPE) approximation, comparing the full (partially degenerate) results to the nondegenerate 
and degenerate limits. The top and bottom panels are for density p = lO^^gcm"^ and 
p = lO^^gcm"^, respectively. The degenerate curves stop at T = Tp/vr. The inset shows the 
crossover region from the degenerate to the nondegenerate regime. 
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Fig. 5. — = 0) [Equation (40)] versus temperature T based on the OPE approximation 
and from chiral NN interactions at N^LO, with a fitting curve to the central value of the 



band. The N LO band is due to different NN potentials employed (Entem & Machleidt 



2003 Epelbaum et al. 2005). In addition, we show results based on the T matrix for NN 



scattering using phase shifts and mixing angles from NN-OnLine (2011 ) for i^iab ^ 350 MeV. 
For temperatures T > 30 MeV, the T matrix results are shown by a thin line, because 
contributions from higher i?iab to = 0) can become important. 
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Fig. 6. — Spin relaxation rate I/to- versus density p for conditions typical of the outflow phase 
in supernova simulations [given by Equation ([2])]. The full (partially degenerate) results are 
compared to the nondegenerate limit for the OPE approximation. In addition, we show the 
spin relaxation rate based on chiral NN interactions at N'^LO and based on the T matrix for 
NN scattering. Note that the density dependence is linear in the nondegenerate limit [see 



Equation (41) 
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Fig. 7. — [Equation (40)] versus oj/T for temperatures T = 5, 10,50MeV. The results 
are based on the OPE approximation. Thick (thin) hues are for the exact u dependence 
(on-shell approximation, see text). 
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Fig. 8. — Relative difference, To- (on-shell)/ro- (exact) — 1, versus density p based on the OPE 
approximation, comparing the exact u dependence to the on-shell approximation (see text) 
for u = T, 3T, and lOT. The comparison is for conditions typical of the outflow phase in 
supernova simulations [given by Equation ([2])]. 
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Fig. 9. — Spin dynamical structure factor S{uj/T) at long wavelength g = as a function of 
energy transfer uj divided by the temperature T for densities lO^^gcm"'^ (upper panel) and 
lO^^gcm^'^ (lower panel). The corresponding temperatures are given by Equation The 
blue bands show our results based on chiral NN interactions at N^LO, while the solid lines 
show results for the OPE approximation. For p = lO^^gcm"'^, we also compare the non- 
degenerate approximation (dashed line) with the full OPE results, while at lO^^gcm"^ the 
results for the nondegenerate approximation are indistinguishable from the full calculation 
on the scale of the figure. The dotted lines show the OPE results of Hannestad Raffelt 



(1998) [HR 1998] as implemented in the Basel code for core-collapse supernova simulations 



(Fischer et al. 2012). 
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